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Abstract 

We perform a Monte Carlo simulation of two-dimensional A^-step interacting 
self-avoiding walks at the 6 point, with lengths up to A = 3200. We compute 
the critical exponents, verifying the Coulomb-gas predictions, the 0-point tem- 
perature Tq = 1.4986(11), and several invariant size ratios. Then, we focus 
on the geometrical features of the walks, computing the instantaneous shape 
ratios, the average asphericity, and the end-to-end distribution function. For 
the latter quantity, we verify in detail the theoretical predictions for its small- 
and large-distance behavior. 



1 Introduction 



Self- avoiding walks (SAWs) have been extensively studied during the years. They have 
a rich mathematical structure and represent one of the simplest models of critical 
behavior. Moreover, they are also relevant for the understanding of the universal 
features of large- molecular- weight macromolecules in solution [T3-|8]. Indeed, the study 
of the asymptotic behavior of long SAWs (and also of walks with different architecture, 
like self-avoiding rings or stars) allows one to obtain predictions for several structural 
and thermo dynamical properties — for instance critical exponents, structure factors, 
osmotic pressure, etc. — of polymers both in dilute and in semidilute good-solvent 
solutions in the limit of infinite degree of polymerization. These predictions are 
in good agreement with the (usually much less precise) experimental results. By 
adding an attractive interaction SAWs also allow us to study the critical transition 
(^-point) between good-solvent and poor-solvent behavior [5HS]. The phase diagram 
of interacting SAWs is well-known. Above the ^-temperature Tg, the large- behavior 
is temperature-independent and in the same universality class as that observed for 
athermal SAWs, which correspond to the case T = oo. On the other hand, for T < Tg 
typical walks are compact: this is the so-called collapsed or poor-solvent regime. At 
the 6 temperature, interacting SAWs show a quite interesting tricritical behavior. The 
two-dimensional case has been extensively studied by Coulomb-gas techniques and 
conformal field theory (CFT). They have provided exact predictions for the critical 
exponents which have been confirmed by several high-precision numerical studies 
[TUH22]- Also the collapsed phase has been discussed in detail [23H2Z], including the 
crossover behavior close to the 6 point |2B] . 

In this paper we perform a Monte Carlo (MC) simulation of two-dimensional inter- 
acting SAWs close to the 6 point, with the purpose of determining some geometrical 
features of the walks. In particular, we shall focus on the instantaneous shape of the 
walks and on the end-to-end distribution function (EEDF). As is well known, walks 
are not instantaneously spherical. Their shape is usually characterized by consider- 
ing combinations of the eigenvalues of the gyration tensor. One of them, the mean 
asphericity, is an essential ingredient in theoretical studies of the hydrodynamic be- 
havior of dilute polymer solutions [29] - [3T] . The EEDF has been extensively studied 
in three-dimensions, because of its theoretical interest [52H11] . Under the mapping of 
SAWs onto the n — t- a model (or theory), it corresponds to the spin-spin corre- 
lation function, which is the basic object of field-theoretical calculations. Moreover, 
it can be also accessed experimentally, by performing neutron-scattering experiments 
on solutions of end- marked polymers [7]. 

In this paper we will obtain a high-precision determination of the EEDF, which 
will be compared with several predictions obtained by using the standard mapping 
onto the n = spin model and with phenomenological expressions that have been 
shown to be quite accurate in three dimensions. As for the shape of the walks, 
we will determine the average asphericity, comparing it with results obtained in the 
good-solvent regime [15| H6] and for noninteracting random walks [47]. 

The paper is organized as follows. In Sec. [2] we define the interacting SAWs and 
some basic quantities. In Sec. Owe use our MC results to determine the 6 temperature 
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and verify the theoretical predictions for the critical exponents. In Sec. |4]we study the 
walk shape and in Sec.[5]the EEDF. Finally, in Sec. |6]we summarize our conclusions. 



2 Model and definitions 

In this paper we consider self-avoiding walks (SAWs) on a two-dimensional square 
lattice. An A^-step SAW w is a set of A^ + 1 lattice sites Uq = 0, coi, . . ., un, such 
that oji and Wj+i are lattice nearest neighbors. For each walk, we define the energy £ 
as follows: 

N-3 N 

g^-Y^Y,c.^, (2.1) 

i=0 j=i+3 

where 

ifk-a;,| = l; ^2.2) 



Essentially, S is the number of nearest-neighbor contacts without considering the 
trivial contacts between subsequent monomers. 

We consider the ensemble of A^-step walks with partition function 

Zr, = J2e-^', (2.3) 

where /3 is the inverse temperature and the sum is extended over all A^-step walks. 
We will study the behavior close to the 6 temperature P0. For the model we consider, 
the best present-day estimates of (3g on the square lattice are reported in Table [H 
We consider three different observables that measure the size of the walk: 

• the square end-to-end distance 

Rl = {con - ojof ; (2.4) 

• the square radius of gyration 

1 ^ / 1 ^ Y 1 ^ 

=^J^2j:(^^-^^f-^ (2-5) 

i=0 ^ k=0 ^ \ ' J jj=o 

• the square monomer distance from an endpoint 

N 



-]^E(^^-^o)^ (2.6) 



p2 

" i=0 



Correspondingly, we define the universal ratios 



_ {Rg)N _ {RJn ^ _ {Rg)N 
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year 


method 


Pe 


Ref. [TO] 


1987 


EE 


0.75 


Ref. dU 


1988 


MC 


0.65(3) 


Ref. [12] 


1990 


EE 


0.67(4) 


Ref. [13] 


1990 


MC 


0.658(4) 


Ref. [H] 


1992 


EE 


0.657(16) 


Ref. [15] 


1993 


MC 


0.658(4) 


rlei. [Id] 


1 CiCiA 




n Ran/ ^\ 


Ref. [n] 


1995 


MC 


0.665(2) 


Ref. [18] 


1996 


MC 


0.664, 0.666 


Ref. [ig 


1997 


MC 


0.667(1) 


Ref. [ig 


2009 


MC 


0.664(8) 


this work 


2010 


MC 


0.6673(5) 



Table 1: Estimates of Pe on the square lattice. EE stands for exact enumeration, MC 
for Monte Carlo. 



and the combination (the exponents 75/ and are defined below) 

2 \ 1 23 1 

An-2Bn + - = —An-2Bn + -. (2.8) 



le + 2uej ^'2 8 2 



= 

For two-dimensional non-interacting SAWs it has been proved [501I5T] that the corre- 
sponding Fn (with the appropriate exponents 7 and u) vanishes in the limit N ^ 00. 
It has been conjectured and verified numerically [IE] that Fqo = also at the 9 point. 
In a neighborhood of P0 the radii have a scaling behavior of the form 

{R')^ = N'^of{x) x^N^i^-Pe), (2.9) 

with /(O) 7^ 0. More precisely, this scaling form is valid for /3 — )• — > c>o at fixed 
X. In two dimensions, CFT and Coulomb-gas techniques allow us to compute the 
universal exponents (p and ug. They are given by [9] 

^e = ^, (2.10) 

The crossover exponent can be measured directly ai (3 = (3e by considering the 
specific heat h^, which scales as 

h^^^{{£')^-{£)y)^N'*-\ (2.11) 

or the temperature dependence of the radii, 
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AT 

i V 








\^)n 


llN 


100 
800 
1600 
3200 


43.438(8) 
472.72(30) 
1050.19(37) 
2338 7fl 5) 


241.87(7) 
2604(3) 
5784(3) 

12895fl3) 


123.02(3) 
1332(1) 
2957(1) 
6584f5) 


46.304(2) 
429.712(26) 
877.717(21) 
1780 259f53) 


68.10(2) 
909.53(60) 
2034.49(68) 
4460 2f2 5) 


N 


An 


Bn 


Cn 


Dn 


Fn 


100 
800 
1600 
3200 


0.179592(61) 
0.18152(22) 
0.18156(12) 
0.18137(21) 


0.50863(18) 
0.51139(64) 
0.51123(35) 
0.51060(63) 


0.35309(97) 
0.35496(35) 
0.35514(19) 
0.35521(35) 


2.2793(15) 
6.826(15) 

9.54(1) 
13.170(29) 


-0.0009(4) 
-0.0009(14) 
-0.0005(8) 
0.0002(14) 



Table 2: Estimates of {Rg)N, {RI)n, {Rm)N, {^)n, hN, An, Bn, Cn, Dn, and Fn- 



In the analysis of the end-to-end distribution function we will also consider the ex- 
ponent 761, which controls the large- iV behavior of the partition function at the 6 
temperature, 

Zn ~ f^^N^'-^ (2.13) 

where /x is a lattice- and model-dependent constant. The exponent ■jg is universal; 
CFT and Coulomb-gas calculations predict [9] 7e = 8/7. 

3 Determination of the critical exponents 

We performed a MC simulation using the extended reptation algorithm discussed in 
detail in Refs. [52l|53]. We fixed (3 = 0.665, which is the estimate of /3g presented in 
Ref. ^7\, and performed runs for = 100, 800, 1600, 3200. Some results are reported 
in Table |2l Since Pg is not exactly known, we also computed several quantities for 
P = 0.665 — 0.0005n, n = 1,2, 3, 4, using the standard reweighting method. 

In order to determine ug we perform fits of the radii. If we fit all data at Pg = 0.665 
to aN'^'^" we obtain ug ^ 0.573, while, if we discard the results corresponding to 
= 100, we obtain ug ^ 0.577. These estimates are close to the theoretical prediction 
ug = 4/7 ~ 0.5714.... The slight discrepancy is probably due to the fact that 
P = 0.665 is slightly smaller than the 6'-value Pg, so that we are seeing the beginning 
of the crossover towards the good-solvent value u = 3/4. A better analysis consists 
in fitting the data at P = 0.665 and the reweighted data to Eq. ( 12. 9p . Since x is small 
for our data, assuming that the function f{x) is regular at x = 0, we can expand it 
in powers of x. At first order we obtain the scaling form 

{R^)n = aN^"' + bN^'o+'/'iP - Pe), (3.1) 

valid for x = N'^{P — Pg) <^ 1. We first perform fits taking a, h, ug, 0, and Pg as free 
parameters. The results are reported in Tabled as a function of A^min, the minimum 
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N ■ 

'mm 


V 





P9 


(3e (f.exp.) 


Rg 


1 nn 

iUU 

800 


U.O ( ZXjyZ, J 

0.5721(2) 


v.^ovy^ ) 
0.480(5) 


0.6669(1) 


U.OO 1 oyij 

0.6675(1) 


Re 


100 
800 


0.5678(3) 
0.5683(3) 


0.480(5) 
0.480(5) 


0.6677(1) 
0.6676(1) 


0.6670(1) 
0.6670(1) 


Rm 


100 
800 


0.5706(3) 
0.5711(3) 


0.478(6) 
0.478(5) 


0.6672(1) 
0.6670(1) 


0.6671(1) 
0.6671(1) 



Table 3: Estimates of the exponents and of /3e obtained by fitting the reweighted 
data to Eq. (13. ip . The critical value (3g in the rightmost column [(3g (f.exp.)] is 
obtained by fixing the theoretical values for the exponents, i.e., by fitting the data to 



length allowed in the fits. The results obtained in the analysis of the three ratios are 
reasonably close and indicate 

ue = 0.570(2), (3.2) 
(j) = 0.479(6), (3.3) 
/3e = 0.6673(5), (3.4) 
Te = l/(3e = 1.4986(11). (3.5) 

These estimates are the average of all fits and the reported error is such to include 
all results and the corresponding errors. The exponent ug is in good agreement with 
the Coulomb-gas prediction = 4/7 ~ 0.571. On the other hand, the exponent (p 
is significantly larger than = 3/7 ~ 0.429. This may be due to neglected scaling 
corrections and/or to the neglected terms in the expansion of the function f{x), 
although, we must admit, we have no evidence of corrections in our results, which are 
stable with respect to A^min and to the observable considered. 

The estimate of (3g is in good agreement with the most recent one obtained in 
Ref . , in which much longer walks were used. The fit in which we assume the the- 
oretical values of the exponents gives estimates of Pg that have a significantly smaller 
statistical error, see Table [31 However, note that we neglect here scaling corrections 
which can — and probably do, given the observed discrepancy for the exponent (j) 
— give rise to systematic deviations which are larger than the tiny statistical errors. 
Hence, we shall keep the conservative estimate (13. 4p . 

In order to obtain independent estimates of the crossover exponent, we also analyze 
Dj^ and h^. By fitting the MC data we obtain = 0.450(4) (from hjy) and (j) = 
0.436(5) (from -Dat). They differ significantly from the estimate 03.31) . and thus provide 
evidence that the apparent stability observed in the fits of the radii and, therefore, 
the relatively small error, should not be trusted. The estimate from the analysis 
of and Dj^j are in better agreement with the theoretical value than the estimate 



6 




0.5 1 1.5 2 2.5 3 3.5 



Figure 1: Rescaled distribution functions {qa) NPa,N{,(la) for the two eigenvalues of 
the gyration tensor. 



(13. 3p . The still present tiny discrepancies indicate that corrections to scaling and/or 
crossover effects are relevant and give rise to systematic deviations, which are larger 
than the statistical errors. This is consistent with the results of Ref. [T7] which 
found (j) = 0.435(6) with significant scaling corrections. Other recent estimates of 
are = 0.419(3) [20j, = 0.436(7) (in a model with exphcit solvent) pi], and 
(j) = 0.422(12) from the analysis of the partition-function zeroes [22] . 

Finally, we consider the invariant ratios A^r, i?jv, Cn, and F^. MC estimates at 
P = 0.665 are reported in Table [21 If we exclude that data with = 100, they are 
constant within error bars, indicating that scaling corrections and crossover effects 
are smaller than the statistical errors. Conservatively, we estimate the asymptotic 
value by averaging the results with > 1600. We obtain 

Aoo = 0.18151(10) (3.6) 

= 0.51106(31) (3.7) 

Coo = 0.35516(17) (3.8) 

= -0.0003(7). (3.9) 

The ratios Aoo and B^o are in good agreement with the results of Ref. [16]: A^o = 
0.180(1), Boo = 0.510(2). The estimate of F^o is fully compatible with zero, providing 
further support to the conjecture F^o = of Ref. [16] . 



4 Gyration tensor and asphericity 

It has been known for many years that polymers are not instantaneously spherical in 
shape [29] . In order to characterize the shape we consider the gyration tensor defined 
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{(1i)n 


{Q2)n 




Sl 


A 


100 


35.543(8) 


7.8954(8) 


0. 


3182(2) 


0.37668(8) 


800 


386.85(30) 


85.871(34) 


0. 


3183(8) 


0.37569(28) 


1600 


859.93(37) 


190.262(41) 


0. 


3188(5) 


0.37653(16) 


3200 


1916.6(1.5) 


422.06(17) 


0. 


3195(8) 


0.37805(28) 



Table 4: MC estimates of the eigenvalues of the gyration tensor, of the shape factor 
Sl, and of the mean asphericity. 



by 

1 ^ 

QN,al3 = 2('A^ -L 1)2 ~ ^i,a)(^i,/3 " ^j,!^), (4-1) 

which is such that TtQn = R"^. The tensor QN,ai3 is symmetric and positive definite, 
hence it has two positive eigenvalues qi > q2- For — > oo they are expected to scale 
as 

(ga)jv ~ 5,iV2-^ (4.2) 
and to obey a scaling law of the form 

PaAQa) = J^F^ (-^) . (4.3) 
[QajN \{qa}Nj 

In Fig. [1] we report the rescaled eigenvalue distributions. All points fall quite nicely 
onto two different universal curves, confirming the validity of the scaling form ( 14. 3p . 
Note that, although the average eigenvalues differ approximately by a factor of 4.5, 
see Table HI the two distribution functions are similar. They have a sharp peak for 
qa/{<la)N ~ 0.55 (a = 1) and 0.75 (a = 2), a long tail, and go to zero sharply for 
qa/{qa)N ~ 0.25. 

To characterize quantitatively the shape of 6 walks we introduce the shape factors 
^1 = 7^i^' ^2 = = 1 - Sl, ri2 = — — = , 4.4 

{Rg}N {RfiN {q2}N 1 - Sl 

and the mean asphericity 

where q = Xla^'a/^ = For a disk we have Si = 1/2, ^ = 0, while a rod gives 

Sl = 1, ^ = 1. 

At variance with the ratios A^v, • • •, the quantities Si and A aX f3 = 0.665 show 
a systematic drift with A^, which may be an indication of the crossover towards the 
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asymptotic good-solvent value. To take it into account, we use the expected scaling 
behavior close to the 9 point: 

s,,A = m-Pe)N% (4.6) 

Expanding the function f{x) to first order we obtain 

Si,A = a^ + a2{P- Pe)Nt (4.7) 

This implies that we should fit our data to ai + hN"^ . The parameter ai gives the 
0-point estimate of the universal ratio. Using = 3/7, we obtain si = 0.8179(4) if 
we fit all data, and Si = 0.8169(20) if we discard = 100. Conservatively, we quote 
as final result 

si = 0.817(2). (4.8) 
From Si we can compute the ratio of the average eigenvalues: 

ri2 = 4.46(6) (4.9) 

Note that the walks are elliptical, with the major axis being approximately a factor 
of two longer than the minor one. 

The N dependence of the average asphericity is not monotonic, and thus the 
fitting form (14. 7p cannot describe the data up to = 100. We thus only fit the data 
satisfying > 800, obtaining 

A = 0.3726(7). (4.10) 

Note that the error is purely statistical and thus it does not include the systematic 
uncertainty due to the scaling corrections. It is interesting to compare the results for 
the asphericity with those obtained under good-solvent conditions [1^H5]: 

^GS = 0.503(1). (4.11) 

Clearly, at the 6'-point walks are more symmetric than in the good-solvent regime. 
Note that our result is closer to the random-walk value [17] 

^RW = ^ - ^C(3) ^ 0.3964. (4.12) 

Thus, also in two dimensions, 9 point interacting SAWs can be reasonably described 
by random walks (in three dimensions interacting SAWs are effectively random walks 
in the limit A^ — )■ oo, although for finite A^ there are quite strong logarithmic correc- 
tions), as far as the shape is concerned. 



5 End-to-end distribution function 
5.1 Definitions 

If CAr(r) is the number of SAWs starting at the origin and ending in r, we define the 
normalized end-to-end distribution function (EEDF) as 

PNir) = . (5.1) 



9 



The mean squared end-to-end distance is related to -Pjv(r) by 

(i?e)iv = E|r|'^^W- (5.2) 

r 

In most of the studies of the EEDF one usually defines a correlation length ^ which 
is trivially related to Re- 

e = 1{RI)n. (5.3) 

In the following we shall always use ^ to characterize the polymer size. In the limit 
N — )■ oo, |r| — oo, with IrjiV^'^ fixed, the function -P/v(r) has the scaling form 

PN{r)^^J{p)[l + OiN-^)], (5.4) 
where p = r p = \p\, and A is a correction-to-scaling exponent. By definition 

K> 

271 pdpfip) = 1, (5.5) 
2'Kp'dpf{p) = 4. (5.6) 



Several facts are known about /(p). For large values of p it behaves as 

/(p) ^/ooP'^exp(-Dp^) , (5.7) 



where a and 5 are given by 



5 = -^ = 1^2.33, (5.8) 

1 — V0 6 

a = 1^:^ = 0. (5.9) 

2 1 -ve) ^ ' 



For p — > 0, we have [35l|37] instead 



/(p) ^ /op', (5.10) 



where 



<)=2l^ = i. (5.11) 

V0 4 

For the purpose of computing D and 5 from Monte Carlo data, it is much easier to 
consider the "wall-wall" distribution function 

P^,N{^) = Y,PN{x,y), (5.12) 

y 

which represents the probability that the endpoint of the walk lies on a plane at a 
distance x from the origin of the walk. In the large- iV limit, Pw,n{x) has the scaling 
form 

P^A^) = ^fM [1 + OiN-^)] p = M (5.13) 
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Figure 2: The wall- wall EEDF: rescaled combination CPw,Nix) vs. p for different 
values of A^. 



For large p we have 



where 6 is given by (15.81) . D is the same constant appearing in Eq. (15. 7p . and [S] 



5.2 Monte Carlo study 

We studied the EEDF following closely the strategy employed in Ref. to analyze 
the same quantity for three-dimensional non-interacting SAWs. Since our runs were 
performed at /3 = 0.665 < f3g, in principle we should reweight the MC data to 
obtain the EEDF at the 6'-point. This correction is apparently negligible compared 
to the statistical errors, so that we directly analyze the results for /3 = 0.665 without 
additional corrections. 

First, we consider the wall- wall distribution Pw,n{x). In Fig. [2] we report the 
scaling combination C,Pw,n{x) versus the scaling variable p. The scaling is essentially 
perfect on the scale of the figure, confirming the correctness of Eq. (15.131) . Then, 
we study the large-p behavior with the purpose of verifying the asymptotic behavior 
(I5.14P and of determining the constants D and /^,oo- Our data are not precise enough 
to determine and thus we always fix its value in the numerical analysis. We fit 
the data to 




(5.14) 




(5.15) 




In ^P^^N 



ln/^,oo - Dp\ 
ln/^,oo - Dp^- 



(5.16) 
(5.17) 
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N 


Pmax 


Pmin 


6 


D 


100 


4 


2 


2.517(3) 


0.127(1) 






3 


2.501(21) 


0.131(5) 




5 


2 


2.466(6) 


0.142(1) 






3 


2.504(13) 


0.133(3) 




6 


2 


2.565(14) 


0.117(3) 






3 


2.601(29) 


0.109(6) 


800 


4 


2 


2.394(12) 


0.148(3) 






3 


2.410(82) 


0.144(20) 




5 


2 


2.318(9) 


0.173(3) 






3 


2.305(29) 


0.177(9) 




6 


2 


2.378(20) 


0.153(6) 






3 


2.361(43) 


0.158(13) 


1600 


4 


2 


2.426(5) 


0.141(1) 






3 


2.390(30) 


0.150(8) 




5 


2 


2.353(7) 


0.164(2) 






3 


2.358(21) 


0.162(6) 




6 


2 


2.414(15) 


0.144(4) 






3 


2.430(32) 


0.139(9) 


Qonn 
ozUU 


4 


Z 


ZADOyii ) 


u.iooyZ j 






3 


2.423(77) 


0.142(18) 




5 


2 


2.371(15) 


0.159(4) 






3 


2.375(45) 


0.158(13) 




6 


2 


2.469(31) 


0.132(8) 






3 


2.493(67) 


0.125(16) 



Table 5: Estimates of D and 6, obtained by fitting the wall-wall EEDF to Eq. f l5.16p . 



The two fits give correct estimates of D and 6, but only the second one provides a 
correct estimate of fw,oo- We consider only the data belonging to the range pmin < P < 
Pmax- An upper cut-off pmax IS uecdcd since scaling corrections and numerical errors 
increase as p increases. First, since Pw,n{x) = for |x| > A^, the deviations from 
fl5.16p and fl5.17p at fixed become infinitely large as p — )■ N/^ ~ N^~'^'' . Second, 
since the EEDF decreases rapidly with p, for p large there is very limited statistics, so 
that Pw,n{x) has a very large error. But large-p data dominate in the fits, providing 
completely unreliable estimates of the fit parameters. 

Results obtained from fits to Ansatz fl5.16p are reported in Table as a function 
of N = 100, 800, 1600, 3200. The results for 6 do not show systematic dependences 
on the fit parameters pmin and pmax (at least in the range we consider), while they 
show a tiny dependence on A^: apparently, for A^ > 800, S increases with increasing 
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Pmin 


6 


D 


log fw, oo 


100 


2 


2.527(4) 


0.125(1) 


-0.640(5) 




3 


2.545(13) 


0.121(3) 


-0.670(2) 


oUU 


9 
Z 






— u.o ( 1 yi-o) 




3 


2.345(30) 


0.162(9) 


-0.520(53) 


1600 


2 


2.412(7) 


0.145(2) 


-0.607(9) 




3 


2.397(22) 


0.149(6) 


-0.581(37) 


3200 


2 


2.431(15) 


0.140(4) 


-0.620(18) 




3 


2.417(46) 


0.144(12) 


-0.598(78) 



Table 6: Estimates of D and 6, obtained by fitting the wall-wall EEDF to Eq. f l5.17p 
with Pmax = 5. 



N. The reason is not fully clear but it may be again an effect of the crossover towards 
the good-solvent value 5 = 1/(1 — z/) = 4. In any case, the results are reasonably 
consistent with the theoretical prediction 6 = 7/3 ~ 2.333. The constant D varies 
roughly between 0.13 and 0.17 for N > 800, so that we can estimate D = 0.15(2). 

In order to estimate fw,oo, we cannot neglect the multiplicative factor p'^™ and 
thus only the results of the second fit are relevant. From the data at iV = 1600 we 
obtain 

log /^,oo = -0.60(5), /^,oo = 0.55(3), (5.18) 

where the error takes into account the estimates obtained by using all values of A^. 

Let us consider now the radial distribution P/v(r). Such a quantity is not well 
suited for a numerical determination of the scaling function /(p), because of fluctu- 
ations due to the lattice structure. In order to average them out, we will employ a 
procedure already used in this context in Refs. [iOlHTi m]. 

We shall consider two different averages 

^'"'''■■-) = S ^'"('•). (5.19) 



-2 i<r2<r2^ 



Here ri„ = tq + nA and '"in ~ + where ro and A are fixed parametersPj and 
^i,n('"i,n) and N2^ni^'2,n) the number of lattice points with the same paritja of 
that lie in the considered shell. For practical purposes, we measure A in units of the 



^ Procedure ()5.19|) corresponds to fixing the width of the annuh on which the average is computed, 
while (|5.20p corresponds to fixing the area (hence the number of points). 

point {x,y) is odd (resp. even) if a; + y is odd (resp. even). 
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Figure 3: Rescaled EEDF ^^Pi"^^^ with A = 1/15 vs. p. 



correlation length: we define A = A/C, and keep it fixed for all values of N. For A 
fixed, -P/v(r), p|^''(r), and p2^j^{r) have the same scaling behavior as iV — )■ oo. The 
same holds for fixed A, as long as A -C 1. In Fig. |3]we report the rescaled EEDF 
obtained by using the average f l5.19p with A = 1/15. All points fall on top of each 
other, except those with = 100 (this is particularly evident for p < 1.5). This 
confirms the validity of the scaling relation (15. 4p . 

Let us now again consider the large-p behavior. In order to determine the param- 
eters, we perform fits of the form 

log[e'p(^-)] = log/oo- V (5.21) 

for each N and for several Pmin, Pmax, and A. Note that in this case theory predicts 
a = and thus this fit allows us to determine /oo, too. 

Results for > 800, A = 1/15 and the fixed- width average (I5.19P are reported in 
Table [3 The results for 6 are fully consistent with the theoretical value, while those 
for D give roughly D = 0.16(2), which is in agreement with the estimate obtained 
by using the wall- wall EEDF. Estimates using A = 1/5, or obtained by using the 
average (I5.20p give similar results. As for /oo we estimate log/oo = —2.60(15) and 
/oo = 0.082(11). 

More precise estimates of D, /oo, and fw^oo are obtained by fixing S to its theoretical 
value 5 = 7/3. From the analysis of the wall- wall EEDF, using the fit function (I5.17p . 
we obtain 

D = 0.1668(3) /^,oo = 0.625(4), (5.22) 
while from the radial distribution function we have 

D = 0.1656(3) /oo = 0.088(2). (5.23) 
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Pmin 


6 


D 


log/oo 


800 


2 


2.355(32) 


0.158(10) 


-2.566(75) 




2.5 


2.340(42) 


0.163(14) 


-2.52(11) 




3 


2.323(56) 


0.170(20) 


-2.46(17) 




3.5 


2.296(83) 


0.180(32) 


-2.35(30) 


1600 


2 


2.381(13) 


0.150(4) 


-2.610(32) 




2.5 


2.375(17) 


0.153(5) 


-2.588(47) 




Q 
O 




v.Looy I ) 


— A.OOoy iv) 




3.5 


2.351(32) 


0.161(11) 


-2.50(11) 


3200 


2 


2.368(35) 


0.154(11) 


-2.589(80) 




2.5 


2.355(47) 


0.159(15) 


-2.55(12) 




3 


2.340(63) 


0.164(21) 


-2.50(19) 




3.5 


2.317(95) 


0.172(34) 


-2.42(32) 



Table 7: Fit results for the large-distance behavior of the radial EEDF obtained with 
the fixed-width average (I5.19p . Here pmax = 7, A = 1/15. 



The estimates of D obtained in the two cases differ by two combined error bars, 
indicating that the errors are underestimated by a factor of at least two. Multiplying 
all errors by two, we end up with the final estimates 

D = 0.1662(6) /^,oo = 0.625(8), /oo = 0.088(4). (5.24) 

We finally consider the behavior for p — > 0, performing fits of the form 

log/(p) = log/o + ^logp, (5.25) 

see Eq. f l5.10p . Since Eq. fl5.25p is valid only for p — > and for r — > oo (scaling limit) 
data must be analyzed in a window pmin < P < Pmax- We find stable results only 
for A^ = 3200. For lower values of A^ lattice effects are very strong and Eq. fl5.25p 
does not describe the low-r data. If we write pmm = ''^minA^, we find stable results 
for ninin > 1, 0.15 < pmax ^ 0.30, and A quite small, A ^ 10~^. If the parameters are 
in this range we obtain 

9 = 0.255(10), (5.26) 
/o = 0.081(2). (5.27) 

The result for 9 is in perfect agreement with the theoretical prediction ^ = 1/4. An 
improved estimate of /o can be obtained by fixing 9 to its theoretical value. We obtain 

/o = 0.0810(5). (5.28) 
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AT 

i V 






A/T 


A/T 


A/T 


100 


1.778(1) 


4.422(4) 


13.89(2) 


52.2(1) 


226.8(6) 


800 


1.815(4) 


4.66(1) 


15.31(7) 


60.8(4) 


281(3) 


1600 


1.819(2) 


4.69(1) 


15.49(4) 


61.8(2) 


288(2) 


3200 


1.818(4) 


4.68(2) 


15.40(7) 


61.2(4) 


283(3) 


-> oo 


1.821(4) 


4.70(3) 


15.5(2) 


62(1) 


290(10) 



Table 8: The non-trivial even moments M2k,N, k < 6, and the corresponding asymp- 
totic values. 



Finally, we computed the moments 



M2k,N = (5.29) 



see Table El We extrapolated the results by performing a fit of the form 

M2k,N = M2fc,oo + aN-^ (5.30) 
where M2k,oo, a, and A are free parameters. The results are reported in Table [HI 

5.3 Phenomenological expressions 

A phenomenological representation for the function /(p) has been proposed by McKen- 
zie and Moore [35] and des Cloizeaux |7]: 

f{p) ^ /ph(p) = Up'^'^ exp {-D,,,p'^-) . (5.31) 

Here 5ph and 6'ph are free parameters, while /ph and Dp^ are fixed by the normalization 
conditions (15. 5p and (15.61) : 



r[(4 + v)/v] 
l4r[(2 + v)/5ph]j ' 

^ 27rr[(2 + V)Aph] ■ ^^"^^^ 

In three dimensions in the good-solvent regime this expression describes the EEDF 
quite accurately, even taking 6 and p equal to their theoretical value [S]. 

In our case, if we use = 6' = 1/4 and 5ph = 5 = 7/3 we obtain for the two 
constants 

Dph = 0.1794, /ph = 0.06931, (5.33) 
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Figure 4: The EEDF against several phenomenological approximations: (a) we set 
^ph = 7/3 and 6'ph = 1/4 and use Eq. f l5.32p to fix the constants; (b) 6ph and 6'ph are 
determined by fitting the data, while the constants are fixed by Eq. (15.321) ; (c) Sp^, 
6'ph, -Dph, and /ph are obtained by fitting the data. 



which are quite close to the exact results. The resulting curve, curve (a) in Fig. HI 
reasonably describes the EEDF in the large- and small-distance region, but under- 
estimates it in the internediate region 0.2 ^ p < 1.4. As an additional check we 
can compute the invariant ratios M2k- Using the phenomenological expression we 
obtain M2k,ph = 1-77, 4.39, 13.9, 53.1, 237 for k = 2,3, 4, 5, 6. They are not very much 
different from the exact results reported in Table [H the differences varying between 
3% for /c = 2 and 18% for k = 6. Note that discrepancies increase as k increases. 
This is due to the fact that these ratios are increasingly sensitive to the large-p be- 
havior, and the phenomenological expression underestimates the EEDF for large p 
since D^^ > D = 0.1662(6). 

In order to obtain a better approximation, we take 6'ph and 5ph as free parameters, 
fixing always Dp^ and /ph by using the normalization conditions (15.321) . We obtain 



Correspondingly Dph = 0.270, /ph = 0.0795. The resulting phenomenological ex- 
pression describes better the EEDF in the relevant region p < 5, see Fig. HI but 
clearly overestimates the EEDF in the large-p region, given that 5ph is smaller than 
6 = 7/3 = 2.333. As a check we again computed the ratios M2k- For k = 2 we obtain 
^4,ph = 1-825 which agrees with the correct value M4 = 1.821(4) and confirms the 
validity of the approximation for p not too large. However, for k > 3 the obtained 
estimates M2fc.ph are larger than the those reported in Tabled For instance we obtain 
-^6,ph = 4.96 and Mg^ph = 17.65, which overestimate the correct results by 5% and 
14%, respectively. 



V ~ 0.282, 



5ph ~ 2. 



04. 



(5.34) 
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We obtain a slightly better approximation if we keep all constants as free param- 
eters, relaxing the normalization conditions f l5.32p . We obtain 



V = 0.277, V = 1.95, Dph = 0.285, /ph = 0.0805. (5.35) 

The corresponding curve is reported in Fig. HJ as graph (c). It cannot be distinguished 
on the scale of the figure from graph (b), obtained by using parameters ( I5.34p . This 
is not unexpected, since the parameters are quite close to each other. For the choice 
f l5.35p of the constants we have 

POO 

/ 2TTpdpf{p) = 1.042, (5.36) 
Jo 

POO 

/ 2np^dpf{p) = 4.469. (5.37) 
Jo 

The violations of the normalization conditions are therefore reasonably small (4% and 
10% in the two cases). 



5.4 Internal-point distribution function 

As a byproduct of our simulations, we also determined an exponent which is related 
to the internal-point distribution function. We consider the probability P/v M(r) that 
um — 1^0 = r, where um is an internal point, i.e. M < N. In the limit N,M ^ oo, 
r — > oo with rN~^ and M/N fixed, we obtain the scaling expression 

P7v,M(r) ^ ^/int(r/e, M/iV), (5.38) 

where = {RI)n/4: as before. The function /int(p, M/A^) is nonanalytic for p — )■ 0: 

/int(p,M/Ar)~p^^-, (5.39) 

where the exponent 6int is independent of M/N. In two dimensions has been 
computed exactly, obtaining 6'int = 5/6 for noninteracting SAWs and = 5/12 at 
the e point [9]. 

The exponent ^int can be determined by measuring the probability P|^^ that the 
endpoint is a nearest neighbor of the walk. Keeping into account that there are A^ 
internal points we obtain 

It should be noted that this expression only takes into account "distant" contacts. 



since the scaling form fl5.38p is valid only in the limit r — oo. To this nonanalytic 
term we should therefore add the contribution of "local" contacts, which is expected 
to be an analytic function of A^. Thus, we obtain the prediction 
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N 


pENN 


100 


0.63307(4) 


800 


0.71579(10) 


1600 


0.73159(6) 


3200 


0.74407(10) 



Table 9: Probability that the endpoint is a nearest neighbor of the walk. 



At the 6 point this gives 

PS'"'-'-+^+j^. (5.42) 
while for noninteracting SAWs we have 

We have computed P^'^^ at the 6 pointEl and fitted the results with a + b/N^ (see 
Table |9]). The estimates of A allow us to obtain an estimate of 6*1111: 

^int = 0.407(11). (5.44) 

This result is in good agreement with the theoretical value 5/12 = 0.4166 . . . More- 
over, we obtain 

pENN ^ 0.7854(14). (5.45) 



6 Conclusions 

In this paper we present a detailed study of some geometrical properties of two- 
dimensional interacting SAWs at the 6 point. For this purpose we have generated 
walks of length up to = 3200 at /3 = 0.665, which is close to the 6 point value 
f3e = 0.6673(5). 

The main results of this investigation are the following: 

(i) We compute the critical exponents and 0. Our estimate of uq, = 0.570(2) is 
in perfect agreement with the Coulomb-gas prediction p] z/^ = 4/7 ~ 0.571. For 
the exponent 0, we find = 0.479(6) from the analysis of the radii, = 0.436(5) 
from the analysis of their temperature dependence, and = 0.450(4) from the 



■^We have also studied P^^^ for noninteracting SAWs in two and three dimensions. In both cases 
the data are well fitted by a + b/N, which allows us to conclude that 6'int > 2(1 — In two 

dimensions this is consistent with the theoretical prediction, while in three dimensions it implies 

^int > 1.40. 



19 



specific heat (errors are purely statistical). The somawhat large differences 
among these estimates indicate that the neglected scaling corrections are im- 
portant. A reasonable final estimate would be = 0.46(3), which takes into 
account all results with their errors. Thus, we also confirm, although with 
limited precision, the theoretical prediction [9] = 3/7 ~ 0.429. 

(ii) We compute several invariant ratios involving the radii i?^, i?^, and and, in 
particular, we verify a conjecture of Ref. [16]. For — )■ oo, the combination 
Fn defined in Eq. (12. 8 p vanishes, as it does for noninteracting SAWs [501 |5T]. 

(iii) We discuss the shape of the walks, determining, in particular, the average as- 
phericity A. We obtain 

>1 = 0.3726(7), (6.1) 

where the error is purely statistical. Walks are typically elliptic, the ratio of 
the two axes being 2.11(2). For comparison, note that for random walks [H] 
A = 0.3964, while under good-solvent conditions [151I1S] A = 0.503(1). 

(iv) We compute the EEDF. We verify the theoretical predictions for its small- and 
large-distance behavior and provide effective approximations valid in the whole 
relevant range r/C, ^5, that is for r/ (R^) ^ ^ 2.5 (for r > 5^ the EEDF is very 
small) . 
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